HBT interferometry with quantum transport of the interfering pair 
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In the late stage of the evolution of a pion system in high-energy heavy-ion collisions when pions 
undergo multiple scatterings, the quantum transport of the interfering pair of identical pions plays an 
important role in determining the characteristics of the Hanbury-Brown-Twiss (HBT) interference. 
We study the quantum transport of the interfering pair using the path-integral method, in which the 
evolution of the bulk matter is described by relativistic hydrodynamics while the paths of the two 
interfering pions by test particles following the fluid positions and velocity fields. We investigate in 
addition the effects of secondary pion sources from particle decays, for nuclear collisions at AGS and 
RHIC energies. We find that quantum transport of the interfering pair leads to HBT radii close to 
those for the chemical freeze-out configuration. Particle decays however lead to HBT radii greater 
than those for the chemical freeze-out configuration. As a consequence, the combined effects give rise 
to HBT radii between those extracted from the chemical freeze-out configuration and the thermal 
freeze-out configuration. Proper quantum treatments of the interfering pairs in HBT calculations 
at the pion multiple scattering stage are important for our understanding of the characteristics of 
HBT interferometry in heavy-ion collisions. 

PACS numbers: 25.75.-q, 25.75.Gz, 25.75.Nq 

I. INTRODUCTION 

Two-pion Hanbury-Brown-Twiss (HBT) interferometry has been widely used in high energy heavy ion collisions to 
study the space-time geometry of particle-emitting sources [H, IH, [1] . As is well known, the HBT interference occurs 
for a chaotic source and is absent for a coherence source Q . It is therefore important to identify the proper "chaotic" 
source to compare theoretical results with experimental measurements. 

Usual theoretical treatments of the HBT problem argue that the collisions in the thermalisation of pions at the late 
stage of the evolution behaves as collisions of classical particles in a chaotic system, and these chaotic collisions will 
lead to a chaotic thermal freeze-out configuration. According to this argument, the thermal freeze-out configuration 
should therefore be the proper "chaotic" source in HBT measurements and the HBT radii are expected to increase 
substantially with increasing collision energies. However, experimental data indicate that there are only relatively 
small changes of HBT radii when the energy increases from AGS, SPS, to RHIC energies [1, 0, d, [l^. El, [T^. ITsj. 

As the proper identification of the HBT "chaotic" pion source is crucial to the understanding of the phenomenon, it 
is necessary to examine the nature of the pion source during the late stage of its evolution. One envisages that as the 
pion system cools down, the kinetic energies of the pions may initially be quite high so that chemical reactions among 
pions can take place, leading to the conversion of pions into kaons (and vice versa). These chemical reactions will 
lower their intensities significantly when the temperature decreases below a limit. The system will undergo chemical 
freeze-out at which the yield ratios of dominant particles are nearly fixed. Before chemical freeze-out, the identities 
of the particles are still in a state of fiux. Therefore, it is reasonable to take the chemical freeze-out configuration to 
be the initial "chaotic" pion source for the purpose of studying the HBT interference. 

The "chaotic" source is however not a static source from which the pair of interfering identical pions are emitted 
and detected. The chaotic source continues to evolve and the pair of pions must follow the evolution from the chemical 
freeze-out configuration to the state of thermal freeze-out. Only at the state of thermal freeze-out can the pair of 
pions be free of interactions and can be considered "emitted" to the detector. How the pair of pions propagate from 
the initial chemical freeze-out point to the thermal freeze-out point will affect the characteristics of HBT interference 
and will be the subject of the present investigation. 

At temperatures below the chemical freeze-out temperature, pions undergo multiple scattering with the medium 
particles. The scattering can be elastic or inelastic. Thermalisation takes place until the temperature reach the 
thermal freeze-out limit, at which the shapes of momentum distributions of final particles are determined [T^.ll5l[l^. 
Thus, the last stage of the evolution of the pion system makes a transition from a chemical-reaction dominating 
epoch before chemical freeze-out to a multiple scattering dominating epoch between chemical freeze-out and thermal 
freeze-out. 

In the usual treatment of the HBT problem, the last stage of multiple scattering is envisaged as random cascades of 
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classical binary collisions and the spatial configuration of the system after these random binary collisions will continue 
to be chaotic in the newly-evolved and expanded configuration. Such a viewpoint may indeed be valid for some other 
considerations of the state of the system, it is nonetheless inadequate for the HBT interference for a pair of identical 
particles. It is important to note that the HBT interference phenomenon is quantum mechanical in nature, arising 
from the interference of pairs of identical pions traveling through different history paths to reach two detectors. The 
propagation of a pair of pions which undergo multiple scattering with the medium at the late stage of the evolution 
should be treated quantum mechanically [13, [iM [3, US HH, HI] . 

Upon carrying out such an investigation for a pair of test particles in the path-integral method, it was previously 
realized that the collisions of the pion with medium particles at the late stage of the evolution lead to the accumulation 
of the phases of the pion wave functions. The real parts of these phase shifts of the pair tend to cancel while the 
imaginary part leads to an absorption in the two-pion interferometry [itI . [isl [Tol . [20| . As a consequence, the effective 
source size is not necessarily that of the newly-evolved spatial configuration as in thermal freeze-out, but may be 
significantly smaller. We would like to demonstrate such an effect explicitly here by assuming a chaotic source at 
chemical freeze-out and by following the quantum transport of the interfering pair, using the path integral method 
for which the pions need not follow straight-line trajectories. For such a study we use relativistic hydrodynamics to 
describe the evolution of the bulk matter while the paths of the two interfering pions are described by test particles 
following the hydrodynamical fluid positions and velocity fields. Alternative descriptions of particle emission using 
different view points have also been presented previously [U, [2^ . 

In following the evolution of pairs of identical pions, there are additional complications. Pions are produced not 
only directly at chemical freeze-out but also indirectly by particle decays as secondary pions during the period from 
chemical freeze-out to thermal freeze-out. It is necessary to include these secondary pions in our HBT calculations. 
The HBT radius for the evolving source will reflect not only the geometry of the source at a configuration close to 
the chemical freeze-out configuration but the effects of these secondary sources. Because the secondary sources are 
different for AGS collisions and RHIC collisions, we shall study the effects of particle decays for collisions at both 
AGS and RHIC energies. As the present investigation is for demonstration purposes, we shall limit our attention on 
spherical geometry. To assist future investigations utilizing quantum transport of the interfering pair, we present a 
detail description of the procedures how this quantum transport is carried out in the path integral method. 

The paper is organized as follows. In Sec. II, we give a brief review of equations of relativistic hydrodynamics 
and discuss the equations of state for sources produced in collisions at AGS and RHIC energies. We also discuss 
the particle production conditions in this section. In Sec. HI, we review the two-pion HBT correlation function in 
quantum probability amplitudes in a path-integral formalism. We examine the HBT radii for the evolving source, 
and investigate the effects of particle decay and multiple scattering on HBT radii both for AGS and RHIC energies. 
Finally, a summary and discussion are presented in Sec. IV. 

II. EVOLUTION OF THE PARTICLE SOURCE 

As a completely quantum mechanical description of the evolution of the full system is not within reach, we shall 
study the evolution of the bulk strongly-interacting matter by relativistic hydrodynamics, which has been quite 
successful in high-energy heavy-ion collisions. We shall then examine the interference of the pair of pions by following 
their trajectories as test particles in the fluid after the chemical freeze-out. 

A. Equations of relativistic hydrodynamics 

In relativistic hydrodynamics, the d yna mics of the strongly-interacting fluid is defined by local conservations of 
energy-momentum and net charges [23.I25I. [26j . In our model calculations we use conservations of energy-momentum, 
net baryon number, and entropy of the system. The continuity equations of these conservations are (24 . [25l . [26j 

d^T^''ix)^0, (1) 
«(^)-0, (2) 

d^j^ix)^0, (3) 

where x is the space-time coordinate of a thermalize fiuid element in the source center-of-mass frame, T^^'^{x) is the 
energy momentum tensor of the element, j^(x) = n\^{x)u^ and j^{x) = s{x)u^ are the four-current-density of baryon 
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and entropy (rib and s are baryon density and entropy density), and — 7(1, v) is the 4- velocity of the fluid element. 
The energy momentum tensor can be expressed as [24, 25), [2y, 

T'^'ix) = [e(x) +p{x)]u''{x)u''{x)-p{x)g'"' , (4) 

where p and e are the pressure and energy density of the fluid element, and g'^'^ is the metric tensor. 

From the local conservations Eqs. (H])-®, one can get the equations of motion for spherical geometry as 0, 

dtE + dr[{E + p)v] = - — {E + p) , (5) 
r 

dtM + dr{Mv+p) = -—M, (6) 
r 

atA^b + 9r(^b«) - -— iVb, (7) 

r 

dtN, + dr{N,v)^-—N,, (8) 
r 

where E = T°°, M = T""', iVb = nb7, = 57. Using the HALE scheme [H,!!!] and Sod's operator splitting method 
[l^l, one can obtain the solutions of Eqs. ©-([S]) [2Q, [2^, [3l|, after knowing the initial conditions and equations 
of state. In our calculations, we assume that the initial systems are distributed uniformly within a sphere with radius 
To and have zero initial velocity. The grid spacing and time step in the calculations are taken as Ax = 0.04ro and 
At = 0.99Ax. 

B. Equation of state 

In the equations of motion ([S])-®, there are e, p, v, n^, and s five unknown functions. In order to obtain the 
solution of the equations of motion, we need an equation of state (EOS), p(e,rtb,s), which gives a relation for p, e, 
nb, and s |2^. 

1. EOS for system at AGS energies 

At AGS energies, we use a mixed perfect gas of hadrons to describe the particle-emitting source. The number 
density, energy density, pressure, and entropy density of the particle species i for a thermalize fluid element can be 
expressed in the fluid element local frame in terms of temperature T{x) and chemical potential fii{x) as 

n. = ^£ f.E^E^-mjdE, (9) 

~ (27r)3 
1 Ang, 



Pi 



3 (2^)3 



f.E'^E^-m^dE, (10) 

>o 

ME'~m^)^/'dE, (11) 



47rg, 
(27r) 



^ rh/.ln/.T(lT/.) 



where 



HlT.mEjE^~m^dE, (12) 



(13) 



exp[{E - fi,)/T] ± 1 ' 
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and rrii are the internal freedom and mass of particle species i, and the sign (+) or (— ) is for fermions or bosons. 
The fluid energy density e, pressure p, and entropy density s are the sum of e^, pi, and Si for all particle species 
considered in the mixed perfect gas, respectively. The baryon density nb is the sum of Ui for all baryon species in the 
gas. 

For a given set of variable {T, fii, fi2, ■ ■ one can obtain the equation of state p{e,nh,s) from Eqs. ([9|l- (fT2|) 
numerically. Additionally, using HALE scheme to solve equations of fluid dynamics we need the sound velocity of the 
fluid, Cg ~ dp/de, which can be expressed for the mixed perfect gas as 

/>oo 

^ E r ^-'^^jE^-ml {E-^,)E'dE. (14) 



2. EOS for system at RHIC energy 



At RHIC energy, the system undergoes a transition from QGP phase to hadronic phase. The initial baryon density 
in the center rapidity region of the collisions is approximately zero because the total number of pious is much greater 
than the total baryon numbers and there is some tendency of baryon transparency. At zero net baryon density, QCD 
lattice results suggest the entropy density of the system as a function of temperature as 25., JiZ, .33|] 



s{T) 




— tanh — — — 

e+1 I AT 



(15) 



where is the entropy density at the transition temperature Tc, AT (0 < AT < O.lTc) is the width of the transition, 
is a parameter associated with the ratio of the degrees of freedom of the QGP phase to the hadronic phase. For 
AT = 0, the equation of state (|15p reduces to the MIT bag equation of state with bag constant B = (G — l)TcSc/2(0+l) 

[HiSlHl. 

For zero net baryon density we only need to solve Eqs. (O, ([6]), and ([8]). The thermodynamical relations among p, 
e, and s in this case are 



sdT + dp^O, 



(16) 



e = Ts-p. (17) 

From these thermodynamical relations and Eq. (|15p . one can obtain the equation of state p(e, s) and sound velocity 
Cg. At the temperature below (1 — AT)Tc, the source is in hadronic phase and the equation of state is taken as the 
mixed perfect gas of hadrons with zero baryon density. 



C. Particle sources 



How the particle number and its chemical composition may change as the fluid evolves can be described by three 
different possibilities, some more realistic than others. They are all included here to show various limiting cases. As 
the equation of state of the fluid depends on particle number and the chemical composition, they also lead to different 
hydrodynamical evolutions. 



1. Partial chemical equilibrium production 



As the temperature of the fluid decreases, the kinetic energy for the collision between particles decreases. As a 
consequence, dominant chemical reactions of converting pions to kaons will lower their intensities and will eventually 
stop when the temperature decreases below a certain limit. For the problem we are considering, the limiting temper- 
ature occurs when the average collision energy at the chemical freeze-out temperature is equal to the threshold for 
the reaction to change two pions to two kaons. The state of the system when this limit of temperature is reached 
can be described as the state of chemical freeze-out, when the equilibrium yield ratios of the dominating components 
(pions and kaons) are frozen. 
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FIG. 1: Chemical potentials of hadrons for CFO production (dotted lines) and PCE production (solid lines) for AGS and RHIC 
energies. 



Prom the passage of chemical freeze-out to thermal freeze-out, the fluid continues to evolve. Its chemical composition 
can continue to change as particles can decay into other particles with various life-times. Conversely, the decay product 
particles can recombine to form parent particles. Some processes with larger cross sections can also equilibrate even 
below T^^. The chemical composition during this stage is however not in a state of complete chemical equilibrium. 
We can call it the state of partial chemical equilibrium (PCE). This is the most realistic case when we include the 
decay and recombination of particles [3, [sH] ■ We follow Hirano and Tsuda [3] to assume that below T'^^ the chemical 
potential of particle species i, fJ-i{T), obeys 

n,(r,M.) _ n,(r-h^A^f) 

where fii = rii + d j^j rij, and dj^i is the fraction of exited-state particle species j that decays to the stable 

particle species i [ij, [S^ . 



2. Chemical freeze-out production 



For our comparison, it is of interest to examine the less realistic cases when the decay of particles is not considered 
and other chemical reactions do not take place after chemical freeze-out (CFO). In this case, the chemical composition 
of the particles are frozen as at the state of chemical freeze-out and all final particles are produced thermally from the 
source configuration at the chemical freeze-out temperature T'^^ < Tc. Below T'^^ all numbers of hadrons are fixed 
and that the particle number densities obey d^{niU>^) = 0. Eq. is changed to 

n,(r,MO ^ n.(T^^Mf) , . 

.(T,{^j) .(rch,{^f})- ^ ^ 

We can obtain the chemical potential iii{T) as a function of temperature below T"^^ from Eq. ^9]) . which ensures 
keeping the number of the hadron i fixed. 

Figs. 1(a) and (b) give the chemical potential ^i{T) for the sources produced at AGS and RHIC energies, and 
CFO and PCE modes of production. For the CFO production, the chemical potentials of all particles are calculated 
from Eq. p9|) . For the PCE production, the chemical potentials of the stable particles are calculated from Eq. (|18p . 
and the chemical potentials of the excited-state particles are obtained from the relations of chemical equilibrium, 
A*A = Mat + AiTT, A'p = and = 3 x 0.88/i^ 14, 35]. 

In our model calculations for the source produced in collisions at AGS energies, the initial temperature and baryon 
chemical potential are taken as 155 MeV and 540 MeV, which correspond to an initial energy density of eg = 0.56 
GeV/fm'^. For simplicity, we only consider iV, tt, and K as stable particles, and A(1232) as exited-state particle for 
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the source with finite baryon density. The chemical freeze-out temperature is taken as O.STq = 124 MeV, and the 
corresponding chemical potential is 565 MeV. They are consistent with the freeze-out values obtained from hadronic 
abundances at CERN/SPS, NBL/AGS, and GSI/SIS [s^. For the source with zero baryon density produced in the 
collisions at RHIC energy, we take the initial temperature as 250 MeV, which corresponds to an initial energy density 
of eo = 4.64 GeV/fm'^. The phase transition temperature Tc and the width of the transition AT are taken as 170 MeV 
and O.OSTc- For the chemical freeze-out temperature IJ], we take the chemical freeze-out temperature as 0.9Tc — 153 
MeV. We consider tt, and 77 as stable particles and take p and oj as excited-state particles after the chemical 
freeze-out. The value of 8 in the equation of state (fT5|) is taken as 5, by consisting the entropy density calculated 
from Eq. with the result of the entropy density calculated for the mixed perfect gas of hadrons at the chemical 
freeze-out temperature. 



3. Chemical equilibrium up to thermal freeze-out 



Finally, we can also examine another less realistic case when chemical equilibrium can be maintained through-out 
the passage up to thermal freeze-out (TFO), with the emission of particles. Implicit in this scenario is the assumption 
that all through the hydrodynamical evolution, the kinetic energy for the collision between particles is always large 
enough to allow chemical reactions to take place to change the chemical composition. Although this may be a good 
picture for classical fluids when the temperature is comparable to the energy needed for chemical reactions, this may 
not be as realistic for our case when the kinetic energy in these temperatures is lower than threshold energy for the 
reactions of two pions converting into two kaons. 

For the TFO particle production, the state of chemical freeze-out coincides with the state of thermal freeze-out at 
T"^. For AGS energies, the chemical potentials of tt and K are zero and the chemical potentials of N and A are the 
same as the initial baryon chemical potential. For RHIC energy, all final particles have zero chemical potential. 



III. TWO-PION INTERFEROMETRY ANALYSIS 
A. Path integral formalism of the HBT correlation function 



As the identities of particles undergo changes during chemical reactions, the pion source at the moment of chemical 
freeze-out can be characterized as chaotic in nature, for the purpose of investigating HBT interferometry. The system 
subsequently expands, its temperature decreases, and its apparent size increases. During the evolution from chemical 
freeze-out to thermal freeze-out, the dynamics is dominated by multiple scattering of the pions. 

It is in the nature of the HBT interference that under the action of multiple scattering after chemical freeze-out, 
the quantum transport of the interfering pair of pions determines their eventual pattern of interferometry when they 
reach the detector. As a consequence, this quantum transport has an effect on the effective source size observed in 
HBT measurements. We suggest previously that the multiple scattering in the transport leads to an apparent source 
that is close to the original chaotic source before the pair quantum transport, supplemented by absorptions along the 
particle paths [13, [3 iS 113 ■ In order to exhibit the importance of this effect, we review the path integral treatment 
of HBT correlation function under quantum transport as discussed earlier in [l9|. 

For definiteness, we consider the HBT interference of a pair of 7r+ particles and study first the motion of one 
of the two identical pions. As the pion propagates from the state of chemical freeze-out to thermal freeze-out, the 
scattering can be described by short-ranged scalar and vector interactions, v^^^^\q — qi), where q is the coordinate 
of the propagating pion and qi is the coordinate of a medium particle. The propagating pion is also subject to a 
collective flow which can be described by a long-range density-dependent mean-field scalar and vector interactions, 
V^^^\q), as in similar cases in the dynamics of the nuclear fluid [s^l- The Lagrangian for the propagating pion is 
given by 

L{q, q) = ( g, q) + L co\{q, q), (20) 
imf (g, q) = ~ [m. + T/if (,)] v/T^ + q • vl^l (q) - (21) 

LUq, 9) = -Koi (9)71"^ + q • V^:l (q) (9), (22) 

Klr^(9) = E-c:r^('?-*). (23) 

i 

From the Lagrangian of the propagating pion, one obtains the pion three-momentum, 

p = 9L/5q = 7[m^ + V^''\q)](i + V^^) (q), (24) 
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the pion Hamiltonian, 

i?-/ = p- q- ^ = + V^''> {q)] + V'>^^'^ (g), (25) 

and the pion mass-sheU condition, 

[pO _ y^iv) (g)]2 _ [p2 _ (g)]2 _ _^ y(s) (^)]2 ^ (26) 

where 7 = 1 / -^Z 1 — and (q) = V^rnf'"' (9) ~^ ^coi"^ ^'^^ ^ pion produced at x with momentum k to propagate 

in the medium to the thermal freeze-out point Xf and be detected at the detecting point Xd with momentum fc, the 
probabihty amphtude according to the path-integral method is [38, 39] 

K{kx kxd) = I Vq e'^i'^^^kx^-^'i) ^ (27) 



where JVq... is the sum over all paths q from a; to Xd, and the action S{kx — > kxd] q) is 

i'Xd,{path q) 

S{nx ^ kxd\ q) ^ / L{q, €[) dt ^ ~ / p{q') ■ dq' . (28) 



We can separate S{kx — > kxd', q) into different contributions, 

S{kx kxd] q) ^ -k ■ (xd - x) + (5mf (kx, kxd] q) + Scoi{kx, kxd] q), (29) 

/'a;^.(path q) -f ^(pa,th q) 

where 5^i{nx, kxd] q) = - [Prai{q') - k] ■ dq' = - / [Pmf (<z') - k] ■ dq' , (30) 



yxd,(path ij) pXfXpath q) 

5co\{nx, kxd; q) ^ ~ Pcoi(g') ■ dq' = - / Pcoi{q') ■ dq' , (31) 



pcoi(g) = f7Klf(9) + v;°i"^(9), 7Klf('?)q + vl:i(9)V (33) 



Because of the additivity of the collision potentials in Eq. ([23| , the phase shift for multiple collision Scoi is a sum of 
the phase shifts for individual collisions, similar to the case of the Glauber wave function in multiple scattering [17], 

(5coi(Ka;, kxd; q) (5coi,i(Ka;, kxd; g), (34) 

i 

where Scoi,i{nx, kxd] q) is obtained from Eqs. (|3ip and (|33p with the potential Vco\{q — qi) in place of the total collision 
potential Vco\{q)- The propagation amplitude is therefore 

K{kx — > kxd) = Jvqexpl—ik ■ [xd — x) + i6-aii{i^x, kxd] q) + i5co\{i^x, kxd] q)}- (35) 

In this sum over all possible path q in the above equation, we can make the approximation that the dominant 
contribution to the path integral comes from the trajectory along the classical path qc for mean-field motion, as other 
trajectories give fluctuating contributions that tend to cancel each other. Then the amplitude is approximately 

K{kx kxd) ~ exp{-i/c • {xd - x) + iS^nfii^x, kxd] qc) + iSco\{KX, kxd] qc)}- (36) 

It is important to point out that the trajectory in the path integral method is quite general and does not need to follow 
a straight line. Our previous consideration for the propagation of an energetic pion along a straight-line trajectory in 
Ref. [l7| is just a special case of the path-integral method. 
The wave amplitude for the pion to propagate from x to Xd is 

kxd) = A{Kx)e">'«''''^K{Kx kxd) , (37) 

where A{kx) is the production amplitude, 4>o{x) is a random production phase for chaotic source. 
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The probability amplitude for the production of two identical pions K2) at {xi, X2) to be detected subsequently 
as fci at Xdi and ^2 at Xd2 is 

kiXdi)'i'i{n2X2 ^2x^2) + (a^i ^ a;2)|, (38) 

where (a:i <-*■ X2) represents the term symmetric to the former by exchanging xi and X2. The probability P{ki, k2) 
for the detection of two pions with momenta (fci, k2) is the absolute square of the sum of the two particle amplitudes 
(Eq. psp ) from all xi and X2 source points. Because of the random and fluctuating phase 4>o{xi) for a chaotic source, 
the absolute square of the sum of the amplitudes becomes the sum of the absolute squares of the amplitudes 1] . One 
obtains 

P(fcl, fe) = P(fcl)P(fc2) [1 + i?(fclfc2) ], (39) 

where P{k) is the probability of detecting one pion of momentum fc, 

P{k) = \A^{kx)K{kX kXd)\^ = ^e-2Xm5,oi(,^x,fcxf;9e)^2(^^) ^ (40) 

X X 

P{ki)P[k2)R{kik2) = ^ A{k,iXi)A{k2X2)K{kiXi — > kiXdl)K{K2X2 k2Xd2) 
X1X2 

xA{kiX2)A{k2Xi)K*{kiX2 kiXdl)K* {k2Xi k2Xd2)- (41) 

We can evaluate the product K{k.iXi — > kiXdi)K*{K2Xi — > k2Xd2)- It is equal to 

K{k^x^^ kiXdi)K*iK2Xi -> k2Xd2) = e-»'^-(---i)+»'=-(----=) 

^ J "^1 J T^q exp{iSniiiKixi,kiXdi;q) - i5^iiK2Xi,k2Xd2;q') 

+ iScoi{KiXi,kiXdi;q) ~ iS*„i{K2Xi, k2Xd2; q')}- (42) 

The real part of the phase difference Scoiii^iXi, kiXdi] q) — 6*^^{k2Xi, k2Xd2] q') is stationary when the trajectories q 
and q' inside the medium coincide and is random and fluctuating when their trajectory in the medium q differs from 
q' . Thus, we have 

J Vq J Vq' expliSmdniXi, kiXdi; q) - iS^f{K2Xi,k2Xd2]q') 

+ iScoliKlXi, kiXdi;q) - iS*^i{K2Xi, k2Xd2]q')} 

~ j Vqexp{iSnx{{niXi,kiXdi;q) ~ iS*^f{K2Xi,k2Xd2;q) 

+ iSco\iKiXi,kiXdi;q) - iS*^i{K2Xi,k2Xd2;q)}- (43) 

Among all trajectories inside the medium, the dominant contribution to the path integral comes from those along 
classical paths, 

J Vq exp{iSnx{{KiXi, kiXdi; q) - i6*„,i{i^2Xi,k2Xd2]q) 
+ i5coi{KiXi,kiXdi;q) - i5la\{n2Xi,k2Xd2;q)} 

exp{i(5,nf(KlXi, fciXdi;(?c) - i5*ra{{l^2Xl,k2Xd2;qc) 

+ i5coi{niXi,kiXdi]qc) - «<5*oi('«2a;i, A:2Xd2; (Jc)}- (44) 

As correlations occur for ki close to k2 these classical paths are close to each other, the real part of the collisional 
phase shifts approximately cancel and only the additive imaginary part due to absorption remain. We then have 

i5co\{niXi,kiXdi;qc) - i5la\{>^2Xi,k2Xd2;qc) ~ ~ImScoiiKiXi,kiXdi;qc) ~^mScoi{K2Xi,k2Xd2;qc)- (45) 

The imaginary part of the multiple scatte ring ph ase shift 2m6co\ikiXi, kiXdi] qc) represents an absorption along the 
trajectory and can be simplified to be [TtI. [iSl Hol. [20| 

/■a;/, {path q^) 

2Im Sco\{KiXi,kiXdi;qc) = / ?^mod(<7') o^abs rf*?', (46) 
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where UaicdiQ') is the medium density at (g'), aabs is the pion absorption cross section, and Xf is the freeze-out 
point. As correlations occur when fci is close to k2 and the absorptions for fci and ^2 from xi along the trajectory 
inside the medium are nearly the same, JmScoiiKiXi, kiXdi; qc) is approximately equal to Im5coi{,K,2Xi, fea^di; qc). We 
can abbreviate these quantities respectively as Im5coi{ki,xi) and X'mScoi{k2, xi) with the other variables implicitly 
understood. 

For the mean- field part of the phase shift, we have 

/■a:;/, (path Qc) 

Smf{KiXi,kiXdi;qc) - S^f{K2Xl,k2Xd2;qc) = / { [Pl mf (g) - fci] - [^2 mf (<?) ^ ^2] } ' rf? 

'J Xi 

= (pm{{kik2,xi;qc) (47) 

Using / d'^xp{x) to substitute J2x ((40|) and (|4T|) for a continuous pion-emitting source with a density p{x), 

we have 



P{k) = / d'^xpix) e-2^"*-'('='^)A2(Kx) 



(48) 



P(fci)P(fc2)i?(fci,fc2) - / d^a;id4a;2p(a;i)p(a;2)e-2^'"^=°'('=i^"i)e-2^'"*-('=^^"^)|$(xia;2;fcifc2)| 



(49) 



i?(fci,fc2) 



^4^gi(fci-fc2)-2;-|-i0mf(*:ifc2,x;qc)-Im5col(fcl,a;)-Im<5col(fc2,a;)^^^^^^^^ x) 



where 



$(:Eia;2;fcifc2) = -^S^Aixi, Kiki)A{x2, K2k2)e''''-^'+''''-''' + {xi ^ X2) 



A{x,Kk) = A{kx) e^'5mf(«^^'=^f;'?c)^ 



(50) 

(51) 
(52) 



Poff(fclfcl 



P(fcl)P(fc2) ^ 



(53) 



and /{kx) — p(x)A'^(kx) is the pion-emitting function of the chaotic source at x with momentum k that evolves 
asymptotically to fc We obtain the simple result that the multiple scattering of the interfering pair leads 

essentially to a factor containing the imaginary part of the phase shift, which is related to the absorption of the 
particles through the medium [itI. [isl. [l9l, |2(| . Thus, even though the source expands to a greater spatial dimension 
during the late stage from chemical freeze-out to thermal freeze-out, the initial chaotic source at chemical freeze-out 
remains an important element in determining the pair correlation function. 

We shall calculate the absorption factor in subsections IIIC and discuss its effect on HBT results in subsection HID. 
The collective flow (mean-field interaction) 'distorts' the initial momentum k into the final detected momentum k. 
The phase associated with the collective expansion of the source, 0mf (fcifc2j x; qc), is is very small p7| . Our calculations 
indicate that its effect on HBT results can be neglected. 



B. Pion-emitting sources 



In our model calculations for the two-pion HBT correlation functions, the identical pions (for example 7r+) are 
assumed to be emitted thermally from the space-time hypersurface at T'^^ for CFO production and from the hyper- 
surface at T^^ for TFO production. For PCE case, the identical pions includes primary pions produced at the chemical 
freeze-out (CFO production) as well as secondary pions arising from decay products of excited-state particles during 
the chemical freeze-out to thermal freeze-out. The four-dimension density of the pion source for PCE production can 
be written as 



p{x) = n^{x)S{T - t"'^) + ^Dj^^nj{x) , 



(54) 
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FIG. 2: Space-time distribution of the pion-emitting source for PCE production at AGS energies. 
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FIG. 3: Space-time distribution of the pion-emitting source for PCE production at RHIC energy. 

where r'^'* is the chemical freeze-out time in local frame and Dj^.^ is the decay rate Tj times its corresponding fraction 
dj^TT- For example, Da-,-^ = Ta x Dp^^^ = Fp x |, D^^^^ — T^^ x 0.88, and D^o^o^^+^- — VrnT^a{'K^iT^ 
TT+TT^) X 1, where Vr is the relative velocity of the two colliding pions. In our calculations, the values of Fa, Tp, and 
are taken as 120 MeV, 145 MeV, and 8 MeV [4l|. The cross section (T(7r°7r° tt+tt") is equal to the absorption 
cross section, crabs(7r+7r^ tt^tt"), and will be discussed in the next subsection. 

Fig. 2 shows the space-time distribution of the pion-emitting source in our calculations for the PCE production at 
AGS energies. Fig. 3 shows the space-time distribution of the pion-emitting source for the PCE particle production 
at RHIC energy. The produced primary pion source points are distributed near the center of the source and the 
decay product secondary pion source points are distributed in the outer region of the source. One can see that the 
space-time geometry of the sources produced at AGS and RHIC energies are different. 




C. Cross sections of pion absorption 



For the PCE case, the pions propagating in the source during the chemical freeze-out to the thermal freeze-out will 
subject to multiple scattering with the medium particles in the source. Based on Glauber multiple scattering theory 
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[40|, the absorption factor owing to multiple scattering in Eqs. 



and dUl) can be written as [13, [H, [H, Hc 



-2 Tin 5co\ (^-fc) 



exp 



2^ CTabs(7ri) rii((jc) dt{qc) 



(55) 



where means the summation for all particles in the source except for the test pion, crabs(7i'*) is the absorption cross 
section of the pion with the medium particle and ni(qc) is the density of the medium particle i along the classical 
propagation path q^. 

For the source produced at RHIC energy, pions are the dominant component of produced particles. The other 
produced particles considered in our model are if, r?, p, and u. Because the threshold of the reaction tttt KK is 
greater than [yfs^) in the source and the probability of three-pion collision for tttttt — > is very small, the dominant 
contributions for the multiple scattering absorption are from reactions tt'^tt — > p and tt+tt^ —^ tt^tt'^. In our model, 
the cross sections for the reactions tt+tt p and tt+tt" 7r° are calculated by [42.. i43i] 



, 47r 2 r 

o'abs(7r"^7r ^ p) = sin 6i 

Pcm 



(56) 



and 



('■abs(7r^7r 



^V) = ^^sin2(5o-<52), 

9 Pcm 



(57) 



where Pcm = \/ s-mr — 4m^/2 and the phase shifts Sq, Si, and 62 are given by Ref. [43[. 

For the source produced at AGS energies, the dominant contribution for the multiple scattering absorption is from 
the process tt^N ^ A. We calculate the absorption cross section in our model by 



<Tabs(7r+A^ ^ A) 



ao (rA/2)2 



3(V^-"lA)2 + (rA/2)2' 



(58) 



where tua = 1232 MeV and cto = 200 mb. 



D. Calculation of P{k) and P(fci,fc2) 



After knowing the hydrodynamical solution, we have the space-time distributions of temperature, chemical potential, 
and velocity of the evolving source. Then, we can calculate P{k) and P{ki,k2) for constructing HBT correlation 
function based on Eqs. (^5]) and (|^ in the following steps: 

Step 1: Select space-time coordinate of the first test pion, xi, in the source randomly with the probability of the 
space-time source density p{x) (at T^^ <T < T^^, see Figs. 2 and 3). 

Step 2: Generate the momentum of the pion in local frame, n'l — {E[,k'i), according to the Bose-Einstein distri- 
butions at the point xi, which is a function of the temperature and chemical potential at 0:1. 

Step 3: Boost k'i by the Lorentz transform with the source fluid velocity at xi. Then, we get the momentum of the 
pion in the center of mass frame of the source, ki = (i?2, ^i). 

Step 4: Determine the new space coordinate of the pion propagating along classical path at next time by the 
production of the pion's velocity (3i = Ki/Ei and the step of time Ai. If the temperature at the new point is higher 
than T'^^ we re-select the pion from step 1 to step 3. 

Step 5: Boost k'i by the Lorentz transform with the source velocity at the new point to get the momenta of the 
pion at the point and calculate Wmsi = "Ylj <^a.hs(.'^j)'nj\/ (/JiAt)^ at this new point. 

Step 6: Repeat steps 4 and 5 until the point a;/ at the temperature T^^ to get the final (detected) momentum of 
the pion, fci, and accumulate Wmsi to get the values of e^^-^'"''™i(^i''^i) based on Eq. ([55)1 . 

Step 7: Repeat step 1 to step 6 to get the corresponding quantities for the second test pion. 

Step 8: According to Eqs. and (1^5]) calculate the weights 



-2 Xm (5coi ) 



(59) 



and 



-2 Xm (5coi (xi ,ki) —2 Xm (Scoi (2:2 ■A;2 ) \_ 



X e 



iki ■x\--\-ik2 ■X2 



E[/E,^/E',/E2 
2 

iki ■X2+ik2 ■X2 



(60) 
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FIG. 4: The two-pion correlation functions for the sources produced at AGS and RHIC energies. 

for P{ki) {i — 1, 2) and P{ki, k2), where E'^j ~ 1, 2) is the energy of the ith pion in the local frame at Xj, which 
is obtained from Ei by a reverse Lorentz transformation with the source velocity at xj . 

Step 9: Repeat step 1 to step 8 and accumulate Wi {i = 1,2) and W12 for each (ki,k2) bin, respectively. Then 
we get the one-pion and two-pion momentum distributions P{k) and P(fci,fc2) for constructing HBT correlation 
fimctions. 

E. HBT results 

The two-particle HBT correlation function C(fci, /C2) is defined as the ratio of the two-pion momentum distribution 
P{ki, k2) to the the product of the single-pion momentum distribution P(fci)P(fc2), 

Using the relative momentum of the two selected pions, q = |ki — k2|, as variable, we construct the two-pion correlation 
function C{q) from P{ki, k2) and P(ki)P(k2) by picking the identical pion pairs emitted from the source and summing 
over ki and k2 for each q bin [13, Sj]- Fig. 4 (a) and (b) show the two-pion correlation functions for sources produced 
at AGS and RHIC energies, respectively. In Fig. 4, the symbols circle and up-triangle symbols are the results for 
the CFO and PCE particle productions. The square symbols present the results for the TFO particle production, 
for this case the HBT calculations of quantum path-integral reduce to a usual HBT calculation for a certain source 
configuration. The curves in Fig. 4 are obtained by fitting the model calculation data of the two-pion correlation 
functions with the parametrized formula 

C{q) = 1 + Ae^9'^' . (62) 

Table 1 gives the fitted HBT radii R and chaotic parameter A in this parametrized formula. In our model calcu- 
lations, the initial source radius, ro, is taken as 6 fm both for the sources produced at AGS and RHIC energies. We 
can see that the HBT radii for the CFO and PCE productions are obviously smaller than those for the TFO particle 
production, especially at RHIC energy. This arises because that the thermal freeze-out is the latest states of the 
expanding sources, which have the largest space geometry of the sources, and at RHIC energy the expanding velocity 
is higher. In order to investigate the effects of the excited-state particle decays and multiple scattering absorption 
on the HBT results in the PCE case, we also examine the HBT correlation functions without the multiple-scattering 
absorptions by letting the absorption factor equal to 1. The corresponding HBT results are listed on right column of 
Table 1. One can see that the HBT radii for the PCE without multiple-scattering absorptions case are larger than 
those for the CFO case because of the excited-state particle decays. The multiple scattering of the identical pions with 
the medium particles in the source leads to an absorption factor to the source density for calculating the two-pion 
HBT correlation functions. It constrains the pions produced near the center of the source and leads to a larger space 
configuration of the pion-emitting source. 

By comparing the HBT radii for the PCE case with and without absorption for the source produced at AGS 
energies. Table 1 indicates that the effect of the large absorption by delta formation increases the HBT radius, as the 
source of pions are located farther out in the expanding system. On the other hand, the HBT radius for the PCE 
case is close to that for the PCE case without absorption, for collisions at RHIC energies. This indicates that for 
RHIC collisions, the effects of absorption is small because the cross sections for tt+tt" — > 7r°7r° and tttt p are small. 
Our results show that at RHIC energy the values of A for the CFO and PCE cases are smaller than unit even for the 
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TABLE I: The HBT fitted results (the unit of source radii R is fm). 







PCE 




CFO 




TFO 


PCE non-abs. 


AGS 


J? = 


4.20 ±0.15 


R = 


3.51 ± 0.11 


7? = 


4.64 ±0.16 


i? = 3.79 ±0.13 




A = 


0.97 ±0.05 


A = 


0.97 ±0.05 


A = 


0.99 ±0.06 


A = 0.98 ±0.05 


RHIC 


R = 


4.22 ±0.23 


R = 


3.74 ±0.15 


R = 


7.28 ± 0.42 


= 4.44 ±0.23 




X = 


0.70 ±0.06 


A = 


0.84 ±0.06 


A = 


0.94 ±0.12 


A = 0.79 ±0.07 



chaotic sources we consider in the paper. Investigating the reason of the A results will be of interest in our future 
work. 



IV. SUMMARY AND DISCUSSION 

Quantum transport of the pair of interfering pion pair has an important influence in the properties of HBT in- 
terferometry in heavy-ion collisions. Staring with a chaotic source at chemical freeze-out, we examine the effects of 
this quantum transport using the path integral method, in which the evolution of the bulk matter is described by 
relativistic hydrodynamics while the paths of the two interfering pions are described by test particles following the 
hydrodynamical fluid positions and velocity fields. By following the trajectories of the interfering pion pair after 
chemical freeze-out up to thermal freeze-out, we find that the quantum transport leads to a cancellation of the real 
part of the phase shifts but the imaginary part leads to an absorption of the propagating pions. The transport gives 
rise to an a HBT radius closer to the radius at chemical freeze-out. We examine in addition the effects of particle decay 
from which secondary pions are produced after chemical freeze-out. As the system continue to expand after chemical 
freeze-out, the spatial dimension of the secondary pions are greater than those of the chemical freeze-out source. 
Secondary pions will lead to a greater HBT radius. The combined effects of of multiple scattering and excited-state 
particle decays lead to HBT radii for the PCE case larger than those for the CFO case, but smaller than those of the 
TFO case. 

The nature of the particle-emitting sources produced in high energy heavy ion collisions at RHIC energy is different 
from that produced in the collisions at AGS energies. At AGS energies, the sources may be described by mixed 
hadronic gas with finite baryon densities. The main stable particles in the sources are nucleons and produced pions, 
and the main excited-state particle may be considered to be A(1232). Whereas at RHIC energies, the sources undergo 
a phase transition from the quark-gluon plasma (QGP) phase to the hadronic gas phase, and with an approximate 
zero baryon density. The main stable particles are produced pions and kaons, and there are many excited-state 
particles, such as p and io in the sources. At AGS energies this effect increases the HBT radius for the PCE case to be 
substantially larger than the HBT radius for the CFO case. Whereas at RHIC energies, the HBT radius for the PCE 
case increases only slightly from the CFO value and is about the same as the PCE case without absorption, because 
the absorption cross sections are small. 

In our calculations we consider only the spherical evolving sources for simplicity. In high-energy heavy-ion collisions, 
the evolution of the system is different in the longitudinal direction (beam direction) and the transverse direction. 
This effect is important for sources produced at RHIC energies. Therefore, it will be of great interest to examine the 
HBT interferometry with quantum transport of the interfering pion pair for more reasonable evolving source models 
and investigate the HBT radii in different directions ("out", "side", and "long" directions) in the future. 
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